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Using a recently developed perturbative approach, which considers Hubbard operators as funda- 
mental excitations, we have performed electronic self-energy and spectral function calculations for 
the t — J model on the square lattice. We have found that the spectral functions along the Fermi 
surface are isotropic, even close to the critical doping where the d-density wave phase takes place. 
Fermi liquid behavior with scattering rate ~ uj and a finite quasiparticle weight Z was obtained. Z 
decreases with decreasing doping taking low values for low doping. Results are compared with other 
ones, analytical and numerical like slave-boson and Lanczos diagonalization finding agreement. We 
discuss our results in the light of recent ARPES experiments in cuprates. 

PACS numbers: 71.10.Fd, 71.27.+a 



I. INTRODUCTION 

Since the discovery of high Tc superconductivit}*i a 
large part of the solid state community accepted that the 
t — J model is fundamental for understanding the physics 
of cuprates^. However, in spite of the great deal of 
work done, important questions about this model are still 
open. In particular, the one-electron spectral function is 
one of the most relevant among them. The recent im- 
provement of ARPES experiments have allowed to take 
access to more refined information about one-electron 
spectral functions and self-energy effects^, both consid- 
ered relevant for understanding the physics of cuprates. 

The main problem for calculating spectral properties 
in the framework of t — J model is the treatment, in 
a controllable way, of the non-double occupancy con- 
straint. There are several analytical and numerical ap- 
proaches used to treat the constrained algebra of the t—J 
model. We will mention some of them. From the ana- 
lytical side can be mentioned: a) Self consistent Born 
approximation^ which is appropriate for the one-hole 
problem, b) Slave fermion (SF)^, even if it is accepted 
that the method is reasonable for low doping, there are 
no many calculations of spectral functions which require 
the evaluation of fluctuations above the mean held level, 
c) Slave boson (SB)^, unlike the SF, seems to be ap- 
propriate for describing the metallic regime. However, 
like SF, the treatment of fluctuations above the mean 
field is not trivial (we further discuss this point below). 
From the numerical side can be mentioned: a) Quan- 
tum Monte-Carlo (QMC), which is suitable for calculat- 
ing spectral functions for one-hole case? while for finite 
doping the sign problem makes the calculation uncon- 
trolled, b) Lanczos diagonalizationi and its finite tem- 
perature version^ which are limited to finite clusters. As 
a consequence, there is no a single method covering all 
situations, therefore it is important to complement ana- 
lytical with numerical methods and viceversa. 



In Ref. iQl we have developed, for J = 0, a pertur- 
bative large- A^ approach for the t—J model based on 
the path integral representation for Hubbard-operators 
(or X-operators) which, in what follows, will be called 
PIH method. This method deals with X-operators as 
fundamental objects and is not based on any decoupling 
scheme; thus, there are no complications as gauge fix- 
ing and bose condensation like in the SB approachiS. 
Recentljiii, the PIH approach was extended to the case 
of finite J . The obtained phase diagram and charge- 
correlations were compared with other calculations based 
on SE^^ and Bayn Kadanoff functional theory (BKF) — 
finding good agreement. 

The aim of the present paper is to present one-particle 
spectral function and self-energy calculations using the 
PIH approach. We show that the method is useful for ex- 
plicit calculation of spectral properties, enabling to sum- 
up, systematically, fluctuations above mean field solution 
given reliable results. 

The paper is organized as follows. In section II, after 
a brief summary of the PIH method, we show the an- 
alytical expressions for the self-energy and the spectral 
function. In section HI, the results are compared with 
available SB ones. In section IV, the results are com- 
pared with Lanczos diagonalization ones for different k's 
on the Brillouin zone (BZ) and different doping levels. In 
section V we present a detailed analysis for self-energies 
and spectral functions for several doping levels and J. 
Conclusions and discussions are given in section VI. 



II. BRIEF SUMMARY OF THE FORMALISM. 

SELF-ENERGY AND SPECTRAL FUNCTION 

CALCULATION 

PIH approach was developed extensively in previous 
pg^pgj,g9jiiji4ji5 g^jjjj jjj ^Yi\s section we list the main use- 
ful steps for explicit calculations of the self-energy and 



spectral functions. 

We associate with the A^-component fermion field fp 
the propagator, connecting two generic components p and 



must be determined minimizing the leading order theory. 
From the completeness condition {Y,p Xf'+X'^^ = iV/2), 
ro is equal to 6/2 where 6 is the hole doping away from 
half-filling. The expression for A is 



Gj?'(k,c.n) 



Ek. 



(1) 



which is 0(1). The original spin index a — ± was ex- 
tended to the index p running from 1 to A^. 

The fermion variable fip is proportional to the 
fermionic X-operator X- ^, fip = {l/\/Nro)X- ^ , and can 
not be associated with the spinons as in SB. In Eq. |^, 
Ek [Ek = —2{tro + A){coskx + cosky) — y) is the elec- 
tronic dispersion in leading order, where t is the hopping 
between nearest-neighbors sites on the square lattice and 
/i the chemical potential. The mean field values tq and A 
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where np is the Fermi function, J the exchange interac- 
tion between nearest-neighbors and N^ is the number of 
sites in the BZ. For a given doping 5] fj, and A must be de- 
termined self-consistently from {I — S) = ^ J^k ''^F{Ek) 
and Eq. 10). 

We associate with the six component boson field 
5X'' = {SR , (5A, r^, r?'. A'', Av), the inverse of the 
propagator, connecting two generic components a and b, 
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The bare boson propagator D^^j is 0{1/N). The first 
component SR of the dX'^ field is related to charge fluc- 
tuations by Xf° = Nroil + dR^), where ^0° is the Hub- 
bard operator associated with the number of holes. SX is 



the fluctuation of the Lagrangian multiplier Ai associated 
with the completeness condition, r^ and A^ correspond, 
respectively, to the amplitude and the phase fluctuations 



of the bond variable A^ = A(l- 



iA^) where i] — x,y. 



J 



The three-leg vertex. 
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represents the interaction between two fermions and one boson. 

The four-leg vertex, A^^ , represents the interaction between two fermions and two bosons. The only elements 
different from zero are: 
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a) Propagators and vertices 
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FIG. 1: a) Summary of the Feynman rules. Solid line repre- 
sents the propagator G^"' (Eq.(l)) for the correlated fermion 
fp. Dashed line represents the 6x6 boson propagator D^"' 
(Eq. (3)) for the 6-component field SX"^. Note that the com- 
ponent (1, 1) of this propagator is directly associated with the 
X"" charge operator. A^^ (Eq. (4)) and A^^ represent the 
interaction between two fermions fp and one and two bosons 
SX"" respectively, b) Hab contributions to the irreducible bo- 
son self-energy, c) Contributions to the electron self-energy 
E(k,cj) through 0(l/iV). 
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Each vertex conserves the momentum and energy and 
they are 0{1). In each diagram there is a minus sign 
for each fermion loop and a topological factor. A brief 
summary of the Feynman rules is given in Fig. ^a). 
As usual in a large- iV approach, any physical quantity 
can be calculated at a given order just by counting the 
powers in 1/A^ of vertices and propagators involved in 
the corresponding diagrams. 

The bare boson propagator D^f^ is renormalized 
in 0{1/N). From the Dyson equation, (Dab)^^ — 

(-^afc )~^ — Hab, the dressed components Dab (double- 
dashed hne in Fig. H^b)) of the boson propagator can be 
found after evaluating the 6x6 boson self-energy matrix 
^ah- Hab niay be evaluated by Feynman rules through 
the diagrams in Fig. ^b). 

In the present summary there is no mention of the 
ghost fields. They were already treated in previous pa- 
pers and the only role they play is to cancel the infinities 
given by the two diagrams shown in Fig. H^b). 

From the iV-extended completeness condition it may 
be seen that the charge operator X°° is 0{N), while the 
operators X^^ are 0(1). This fact will have the physical 
consequence that PIH weakens the effective spin interac- 
tion compared to that one related to the charge degree 
of freedom. 



The component Dnji (component (1,1)) of the 6x6 
boson propagator is related to the charge-charge correla- 
tion function x^ by 



X^iq,Vn)^N 



DRR{q,Vri)- 
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In Ref. 13, O we have pointed out that, in 0(1), the 
charge-charge correlation function shows the presence of 
collective peaks above the particle-hole continuum. 

In what follows self-energies and one-particle spec- 
tral functions are calculated by means of the Feynman 
rules. The Green's function ^ corresponds to the A^- 
infinite propagator which includes no dynamical correc- 
tions; these appear at higher order in the \/N expansion. 
For obtaining spectral densities, the self-energy E is cal- 
culated. Using the Feynman rules, the total self-energy 
in 0(1/A^) is obtained adding the contribution of the two 
diagrams shown in Fig. Q^c). The analytical expression 
for S, for a given channel p, results: 



E(k,iw„) = E(i)(k,i 



E^^^(k,icj„ 



(10) 



where 



S(ink,*c.„) = ^Y^KPPDaM.iv^) 



q^i'n 



X GW(k-q,z(a;„-^.„))Af (11) 



and 



E(2)(k,zc.„) = -l^A^^i?,b(q,^^„). (12) 

The sum over repeated indices a and b is assumed. The 
renormalized boson propagator Dab plays a similar role as 
the phonon propagator when dealing with the electron- 
phonon interaction in simple metals. Therefore, in the 
calculation of S(k,w) through 0(1/A^), enter the band 
structure effects and collective effects associated with the 
charge degrees of freedom (see Eq. ©). Using the spec- 
tral representation for the boson field. Dab: we obtain 



E(i)(k,zc^„) = 
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where 



B''\q,iy) = -2 lim Im[i?ab(q, it^„ -^ly + ij^)]. (15) 



After performing the Matsubara sum and the analyt- 
ical continuation iw„ = u + ii], the imaginary part of S 

is 



where ns is the Bose factor, and the 6-component vector 
ha{k,q,v) is 



ImS(k,w) = -^^ha{k,q,uj - Ek-q) 
q 
^'^''(q, uj - Ek-q) hb{k, q,u;- Ek-q) 
x[nF{~Ek-q)+nB{u) -Ek-q)] (16) 
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It is interesting to show the more compact result for 
the case J = 0: 



^ q 

+ 2 n B^'^ifl^u; ~ Ek-q) + B^\(l,u; - Ek-q)} 
x{nF{-Ek-q) + nB{uj-Ek-q)) , (17) 



where J7 — [Ek-q + 2/i + cj)/2. Using the Kramers- 
Kronig relations can be determined ReS(k, w) from 
Iml](k, Lj) and compute the spectral function A(k, w) — 
— ilmG(k, w) as 
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J/t = 0.3 
5 = 0.2 !; 



k = 0.67(71/2,71/2) 
k = 0.86(71/2,71/2) 
k= 1.20(71/2,71/2) 
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TT (w - Sk - ReS(k, tj))2 + Iml](k, tj)2 



(18) 



The self-energy is calculated using the propagator 
G(k, w) for the /-operators which are proportional to the 
fermionic Hubbard operators and then, they can not be 
related to usual fermions. 



III. COMPARISON WITH SLAVE-BOSON 

While many papers on SB have been published on the 
mean field level there are few calculations including fluc- 
tuations above the mean field which are necessary for 
the estimation of spectral functions. This shows that, in 
spite of the popularity of the SB method, it is not triv- 
ial to implement this kind of calculations. Even if PIH 
results have, at the mean field level, a close connection 
with SB (see Refs. [3, [ill) it is relevant to compare both 
approaches beyond the mean field level. To our knowl- 
edge, there is only one paper where spectral functions 



FIG. 2: Solid, dashed, and doted lines are the spectral 
functions for k = 0.67(7r/2, 7r/2), k = 0.86(7r/2, 7r/2) and 
k = 1.2(7r/2, 7r/2) respectively. The frequency uo is in units of 
t which is considered to be 1. For comparison with SB, in the 
inset we have included the results from Ref. [la . 



have been calculated in the framework of SBi^. In that 
paper, their authors present results for J = 0.3 and dop- 
ing 5 — 0.20 for three different k-points on the BZ. 

In Fig. [21 we present PIH spectral functions A{k,uj). 
The calculation was done for J = 0.3, 5 — 0.20 and 
for k = 0.67(7r/2,7r/2), k = 0.86(7r/2,7r/2) and k = 
1.2(7r/2, 7r/2), which are exactly the same conditions of 
Fig. 3 in Ref. [la- In the inset of Fig. [5|we have included 
SB results for comparison. As can be seen, the spectral 
functions have some similarities with SB. We have ob- 
tained a low energy peak and a pronounced structure at 
large binding energy. The low energy peak is the quasi- 
particle (QP). The other features, at large binding en- 
ergy, are incoherent spectra (IS). For k ~ 0.86(7r/2, 7r/2) 
the QP peak crosses w = 0, where the chemical potential 



is located. In spite of similarities between present results 
and those of Ref. [1^ there are some differences: a) Our 
IS is located at binding energy larger than in SB. For 
instance, for k = 0.86(7r/2,7r/2), IS is at tj ^ — 3i while 
the corresponding one in SB is at cj ^ — 2t. b) Our QP 
peak is less dispersing than in SB. As it is well known, 
self-energy effects depresses Fermi velocities (vp) respect 
to the bare one «'"''=), v*p = Zv)?''^, where Z = (1- |§) 
is the QP weight therefore, it is concluded that our self- 
energy effects are larger than in SB. 

In SB there are three different E'si&: E„, Eq and Ei„c- 
In E„, bosons are condensed, in Sq one boson is con- 
densed and the other fluctuates, and in Ei„c both bosons 
fluctuate. These complications are due to the decoupling 
scheme used in SB, so beyond mean field level it is neces- 
sary to convolute spinous and bosons for reconstructing 
the original X-operators. As in the PIH approach there 
is no any a priori decoupling and we work directly with 
the Hubbard operators, we have only one self-energy S 
given by Eq. (|ll|l . A detailed description of the present 
self-energies for different 5 and J is given in section V. 

After comparing spectral functions with some available 
SB results and, in spite of some similarities, we found 
important differences in the self-energy effects between 
both methods. The existence of only a few SB results 
beyond mean field level may be closely related to the 
decoupling scheme which leads to a gauge field theory, 
making hard the implementation of the approach. We 
hope that PIH be useful and a complement of the SB 
calculations. 
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FIG. 3; (Color online) Comparison between spectral function 
calculations obtained by Lanczos diagonalization for all the 
allowed k-points of the 4x4 lattice and results obtained with 
the PIH method. The parameters are J = 0.3 and 5 = 0.75. 
The k-points are presented between parenthesis in each frame. 
The frequency uj is in units of f. In order to plot PIH and 
Lanczos results on the same figure we have used an arbitrary 
scale for the y-axis. 



analysis of the other k's. 

a) Results for doping S — 0.75 



IV. COMPARISON WITH EXACT 
DIAGONALIZATION 

As pointed out in Sec. II, PIH approach weakens spin 
fluctuations over charge fluctuations. For instance, at 
leading order, while there are collective effects in the 
charge channel, the spin channel exhibits the character- 
istic form of a Pauli paramagnet where the electronic 
band is renormalized by correlationsiiii. So that, at 
0(1 /N), the self-energy does not contain collective ef- 
fects, like magnons in the spin channel. Intuitively one 
may think that the method will be better for large doping 
than for low doping. However, the exact role played by 
charge and magnetic excitations as a function of doping, 
in the t — J model, is one of the key points for under- 
standing the physics of cuprates. Many mechanisms have 
been proposed. Some of them privilege charge while oth- 
ers privilege magnetism. 

In order to test the reliability of our results as a func- 
tion of doping, in this section, we compare qualitatively 
present spectral functions with calculations obtained us- 
ing Lanczos diagonalization. For this purpose we have 
performed Lanczos diagonalizationpiS, on the 4x4 lattice 
for 5 = 0.75, 0.5 and 0.3125, with J = 0.3. 

As an example, we will explicitly compare some 
k-points for several dopings leaving to the reader the 



This doping corresponds to 12 holes in the 4x4 lattice. 
There is good agreement between Lanczos and PIH (Fig. 
13)1 for the six k-points allowed in the 4x4 lattice. For 
instance, for k = (0, tt) both calculations show a QP peak 
at around uj ^ 2t and IS at around lo ^ —At. Lanczos 
diagonalization shows a small peak at w ^ 3i which is not 
seen in PIH. However, PIH shows an asymmetric shape 
of the QP peak which can be interpreted as an indication 
of the additional structure observed in Lanczos. For k = 
(0, 7r/2) both, Lanczos and PIH show a QP peak near 
the Fermi level, and IS at a; ^ —At. The additional peak 
that appears in Lanczos at u) ~ I.5t can be associated 
with the non-symmetrical shape of the QP peak seen in 
PIH. 

b) Results for doping 5 — 0.50 

This doping corresponds to 8 holes in the 4x4 lattice. 
Results are presented in Fig. 0] for Lanczos and PIH. 
With decreasing doping from 0.75 to 0.50, both methods 
show more IS. For instance, for k = (0,0) both calcula- 
tions show two well defined peaks below the Fermi level. 
The peak closer to the chemical potential is the QP, and 
the peak near to ~ — 3i is of incoherent character. Both 
methods also present small IS for lo > {). 

For 5 = 0.5, results for J = were presented in Ref. 
[ij in the context of organic materials where PIH spectral 
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FIG. 4: (Color online) Comparison between spectral function 
calculations obtained by Lanczos diagonalization for all the 
allowed k-points of the 4x4 lattice and results obtained with 
the PIH method. The parameters are J — 0.3 and 5 = 0.5. 
The k-points are presented between parenthesis in each frame. 
The frequency uj is in units of t. 



FIG. 5: (Color online) Comparison between spectral function 
calculations obtained by Lanczos diagonalization for all the 
allowed k-points of the 4x4 lattice and results obtained with 
the PIH method. The parameters are J = 0.3 and S — 0.3125. 
The k-points are presented between parenthesis in each frame. 
The frequency u is in units of t. 



functions were also compared with those obtained using 
Lanczos diagonalization as a function of the nearest 
neighbors Coulomb interaction V. 

c) Results for doping S = 0.3125 

This doping corresponds to 5 holes in the 4x4 lattice. 
The results are presented in Fig. Ofor Lanczos and PIH. 
Lanczos and PIH present the QP peak near the Fermi 
level and stronger IS than for S = 0.75 and S = 0.5. 
The increasing of the IS is consistent with the fact that 
the QP weight is Z ^ 0.5, lower than for the previous 
dopings where for S = 0.75 and S = 0.5 is Z ^ 0.9 and 
Z ~ 0.7 respectively. The QP weight will be discussed in 
more detail in next section. 

For k = (0, 0) both calculations show a QP peak 
below Fermi energy. The peak that appears in PIH at 
Lu ^ — 3t is broader in Lanczos and centered at w ^ —At. 
For k — {it, Tr/2) both calculations show a QP peak 
above Fermi level and IS on the top of this peak. The 
well pronounced structure obtained in PIH at uj ^ ~6t 
seems to be missing in Lanczos. Instead, it shows an 
homogeneously distributed IS below the Fermi level up 
to large binding energies of the order of —7t. In the 
frequency range — 5t < lu < —It both methods show 
IS. For k = (0, tt) both methods present a QP near 
the Fermi level and IS at w ~ — 5i. Between those two 
features it is possible to see IS (in the form of several 
peaks in Lanczos) in both calculations. 

For the three studied dopings both methods show that 
while the QP peak disperses through the Fermi surface 
(FS), the edge of the IS moves in opposite direction. This 



result was obtained previously by Stephan and Horscbi^. 
In Ref. [13 the authors also studied spectral functions for 
the t — J model at moderate doping 5 '^ 0.1 by means 
of exact diagonalization. That paper presents strong 
evidences for a large FS for moderated doping levels. 
This result gives an additional support for our bare band 
Ek — —2{t (5/2 -I- A)(cos fca; -I- cos ky) — /i. 

We conclude that PIH and Lanczos results, for the 
three studied 5 values (which cover a broad range of dop- 
ing) , fairly agree considering the different nature of both 
methods. The above comparison gives some confidence 
to our self-energies. The self-energy has additional infor- 
mation such as relaxation times 1/r, quasiparticle weight 
Z (effective mass increasing) which can not be directly 
obtained from Lanczos diagonalization. 

Decreasing doping , Lanczos and PIH both show band 
narrowing. The narrowing is stronger in PIH than in 
Lanczos. For instance, for k = (tt, tt) (Fig. ^, while 
the Lanczos QP peak is at w ~ 4t, the PIH QP peak is 
at w ~ 2i. For k = (7r,7r) (Fig. [SJ, while Lanczos QP 
peak is at w '^ 2t it is at w '--^ It in PIH. For 5 — 0.75 
(Fig. 121) there is good agreement in the QP and IS energy 
positions for each k. This discussion is important in the 
light of ARPES experiments in cupratesSSiSi which show 
Fermi velocities rather independent of doping (see also 
Ref. I23 for discussions) . In PIH the strong narrowing is 
mainly due to the factor 5/2 in the electronic dispersion 
which strongly weakens the i-term. 
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FIG. 6: Quasiparticle weight Z as a function of doping for 
J = 0. 



V. SELF-ENERGY RENORMALIZATIONS 

In this section, we present a detailed self energy and 
spectral functions calculations from PIH. 

Fig. shows the QP weight Z as a, function of doping 
for J = 0. As it was shown in Ref. [ij, for J = the 
homogeneous Fermi liquid (HFL) remains stable for all 
d. For each doping, we have found that the self-energy 
is very isotropic on the FS (see below), making the QP 
weight rather constant. In Fig. Z ^ when (5 — > 0. 
For small 6, Z ^ lA S, which is very close to the observed 
ARPES behavior in LSCCM^. As Z remains finite for 
(5 > 0, present calculation predicts a Fermi liquid (FL) 
behaviour. Fig. |H|also shows that Z ^ 1 when (5 — > 1 as 
expected for an uncorrelated system. 

Let us discuss the case J — 0.3. Fig. [7| shows S, for 
5 ~ 0.3, for three well separated k-vectors on FS. One 
of the k is chosen in the (ll)-direction of the BZ, other 
in the (lO)-direction and the third in between both. The 
upper and the middle panel of Fig. Qshow the ImE and 
ReS respectively. 

As shown in Fig. [7| PIH predicts a rather isotropic 
self-energy on the FS. On the other hand, for each k, 
I](k, w) is very asymmetric respect to w = which can 
be interpreted as a consequence of the difference between 
addition and removal of a single electron in a correlated 



system. Near w = 0, Iml](ki?,w) 



showing FL 



behavior. On the other hand, ReE(ki^,a;) shows, at uj = 
0, a negative slope which is also a characteristic of a Fermi 
liquid. 

Inset of Fig. [7| shows a plot of — ImI](kF, lo) for w < 
for 'kp in the (ll)-direction. We have used t = 0.4 eV^. 
In the range —200 meV < uj < 0. — ImS](k, a;) does 
not saturate as in Fig. 1 of Ref. ^^ The no satu- 
ration of — ImE(k, Lj), up to an energy scale of the or- 
der of —200 (—300) meV, is well established in cuprates 
and clearly it can not be explained by phonons. In ad- 
dition to this feature. Fig. 1 of Ref. |2y shows the 
presence of an additional energy scale of the order of 
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FIG. 7: —Im'E(kF,u!) (upper panel), _ReE(kF,(^) (middle 
panel) and spectral functions (lower panel) for J = 0.3, 
S = 0.30 and three different points on the FS. Solid line, 
dashed line and dotted line are results for the Fermi point 
in the (1, 1), (1,0) directions and between these two respec- 
tively. The frequency u, ReT, and 7mE are in units of t. 
Inset; — I m'S{'kw , i^) for kr in the (ll)-direction in units of 
t = 0.4eV. 



60 — 70 meV which is associated with the kink observed 
in ARPES^. This small energy scale is not seen in our 
— ImS(k,w). Whether the kink is due to magnetic exci- 
tations or additional degrees of freedom like phonons, is 
still controversial'^''^'*. For instance, Yunoki et alm^, us- 
ing variational MC, found no evidence for the kink in the 
context of the pure t — J model and, on the other hand, 
FLEX calculations for the Hubbard model suggest that 
the interaction between QP and spin fluctuations leads 
to the kink (see Ref. 25 and references therein). As men- 
tioned in previous sections the PIH approach weakens 
spin fluctuations over charge fluctuations, which means 
that 1/7V^ self-energy corrections should be calculated in 
order to study the kink, if originated by magnetic exci- 
tations. 

Finally, as expected from the results shown in the up- 
per and middle panels, the lower panel of Fig. |7| shows 
that, for the three mentioned k-vectors, the correspond- 
ing spectral functions are isotropic. 

As it was shown in Ref. IllL in agreement with previ- 
ous calculations^MLiSi, for J = 0.3, the system presents 
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FIG. 8: Quasiparticle weight Z versus doping for the three It- 
points on the FS discussed in Fig. |7| As discussed in the text 
there is no indication of the proximity to the DDW phase. 
The calculations run from (5 > 5c ~ 0.14 where the HFL is 
stable. 



FIG. 9: Q^F(tj) as extracted from Eq. JTZJ for E. As dis- 
cussed in the text, ct'Fiui) is mainly dominated by collective 
excitations of the charge character. The frequency ld and o?F 
are in units of t. 



a flux phase (FP) (also called rf-density wave (DDW)^^) 
for 8 < Sc ~ 0.14. FP was interpreted as a candidate 
for the pseudogap phase of cuprates^SiSS. Thus, it is 
important to study self-energy corrections approaching 
the FP instability from the HFL phase {6 -^ Sc from 
above). Similar calculations to those for 6 — 0.3 show, 
for S > 5c ^ 0.14, isotropic self-energy effects on the FS. 

Fig. [Slshows Z versus S for three k-points chosen as in 
Fig. [7| each one of them on its corresponding FS for each 
doping. Results are for S < 0.3. The QP weight results 
very isotropic on the FS even for doping near Sc- Accord- 
ing to our results the anisotropy, between AT-point ((10)- 
direction) and nodal point ((ll)-direction), observed in 
ARPES spectra in cuprates, can not be interpreted as 
originated by self-energy effects. This is close to the re- 
cent interpretation by Kaminski et al?^ where the scat- 
tering rate was found to be composed by an isotropic in- 
elastic term and a highly anisotropic elastic term which 
correlates with the anisotropy of the pseudogap. In our 
case, the ImE can be interpreted as the inelastic con- 
tribution to the scattering rate and the opening of the 
flux phase, which is mainly of static character-^'-^, as 
the elastic term. (A close comparison with ARPES ex- 
periments in cuprates, which needs a better FS as given 
by the tt' — J model, is in progress). 

It is important to discuss the reason for the non strong 
influence of the flux instability on the self-energy. As 
discussed in Refs. [Til, |2^ flux phase is mainly of static 
and d-wave symmetry character and it is weakly coupled 
to the charge sector. Being our self-energy dominated 
by charge fluctuations (see below), S does not stro ngly 
prove the proximity to the DDW. In terms of Ref. |2^ 
the proximity to the DDW is hidden for the one-particle 
spectral densities. This is in contrast with the self-energy 
behaviour in the proximity of the usual charge density 
wave (CDW) instability. In Ref. [ij it was shown that 
the QP weight Z is strongly affected when the system 
approaches the CDW phase. 



A discussion about the excitations that, interacting 
with electrons, cause the self-energy renormalizations is 
necessary. In the usual many body language, the self- 
energy can be expressed in terms of the relevant quan- 
tity a^ F (bj^ir^, where the notation is chosen in a way such 
that F{u)) gives information on the density of states of 
a boson interacting with the electrons, and a^ about the 
coupling. In usual metals, a^F(uj) contains information 
of the electron-phonon interaction averaged over the FS. 
For simplicity, we will discuss the J = case. Eq. (|17|) 
is conveniently written for reading a^F{i^). In the flrst 
term of the second hand side of Eq. (|17|) we can interpret 
j^BR _ _2ImDi^fl; as the spectral function of the boson 
mediating the interaction, while the remaining squared 
factor, n?, as the coupling. As discussed in Sec. II, 
Drr corresponds to the charge-charge correlation func- 
tion (Eq. (9)). 

Fig. |51 shows a'^F{uj) obtained from the first term 
of Eq. (17). Following the discussion above, a^Fiuj) 
is proportional to the average on the FS of the charge 
densities. Clearly, charge densities survives up to high 
energy causing the large self-energy effects at large uj. 
Since charge densities, in 0(1), present collective peaks 
at the top of the particle- hole continuum^ ^ both, the 
collective excitations and the continuum, contribute 
to a^F{Lo). For instance, the pronounced structure 
at cj ^ ±1.7 in Fig. |51 is mainly due to collective 
fluctuations. The interpretation of the last two terms of 
the second hand side of Ea. l|17|) is less direct and they 
are proper of our strong coupling perturbativc approach. 
However, they are also dominated by collective excita- 
tion of charge character arising from the inversion of the 
matrix D. 

Sum rule: Before closing we will discuss the spectral 
function sum rule J duj A(k, uj). In the framework of the 



t — J model the sum rule is < AT^" 
(1 + (5)/2. Using the relation X^'^ = 



> 



< X° 



>= 



y/Nro fp, in the 



limit N — 2, our sum rule is S (< X"" >), therefore, PIH 
misses a contribution < X°'°' >= (1 — 6)/2 making the 
situation better for large than for low doping. It is impor- 
tant to discuss about a possible origin of this discrepancy. 
As was pointed out in previous papers^SiS, in order to 
guarantee the commutation rules for X-operators not all 
the multiplication rules can be satisfied. For instance, in 
Ref. 32 we have studied the spin system using the four 
X-operators X'^'^ and showed that the formulation leads 
to the well known coherent state path integral represen- 
tation for spins^. The fact that this representation is 
better for lar ge t han for small S (Ref. l34|) was under- 
stood, in Ref. 132 . as a consequence that not all the mul- 
tiplication rules are fulfilled. It is worthy to note that the 
formalism in Ref. l32l reproduces the spinless fermion case 
when they are written using X-operator representation. 
In the present case we deal with the t — J model and in 
order to satisfy the commutation rules, the formalism re- 
quires the constraint X'^^X^''' = x^^X""'^-^'^ which re- 
produces the exact multiplication rule X'^^X'^'^ — X"'°' 
in the limit X°° ^-1 (6^1) making the representation 
better for large than for low doping. In mathematical 
terms, our expansion seems to be appropriate for both 
large N and large 6. This is closely related to the fact 
that the formulation weakens spin over charge fluctua- 
tions. The band narrowing, discussed in section IV, is 
possibly connected with this discussion if a spin term 
< X"''^ > also contributes to the bandwidth. The solu- 
tion of this very hard theoretical problem, and the knowl- 
edge of how important it is as a function of doping on 
different physical quantities, requires not only to make 
an effort on formal level, but also, at the same time, con- 
fronting results with different methods. 



VI. DISCUSSIONS AND CONCLUSIONS 

The recently developed path-integral large- A^ approach 
for Hubbard operators, PIH, was used for calculating self- 
energy corrections and spectral functions, including fluc- 
tuations above mean field solution of the t — J model. 

Similarities and differences with SB were discussed in 
section HI. To gain confidence on our calculation, com- 
parisons of spectral functions with Lanczos results for 



J = 0.3 and for doping 5 = 0.75, S = 0.5 and S = 0.3125 
were performed in section IV. We found fair agreement 
for each k on the BZ. PIH self-energies and spectral func- 
tions for different J and S have been investigated in Sec. 
V. The general characteristics of the self-energy are: 

a) Around tj = 0, Iml](tj) ^ w^ which is characteristic 
of a FL behavior. This is in agreement with the negative 
slope of ReS(w) at w = 0. 

b) ^(iij) is very asymmetric with respect to w = 0, 
indicating the difference between addition and removal 
of one electron in a strongly correlated system. 

c) E(w) has large structures at large negative u; of the 
order of a few t. 

For J = 0, we have shown that Z decreases monoton- 
ically a,s 6 ^ remaining finite for (5 > 0. For small S, 
Z ^ 1.4 5. As expected, in the uncorrelated limit, Z ^ \ 
for (5 — > 1. 

We have also studied spectral functions along the FS 
for different 5 and J = 0.3. For this case, the HFL is 
unstable against a DDW phase for doping 5 < 5c ^ 0.14. 
Since DDW phase was interpreted as a candidate for de- 
scribing the pseudogap state in cuprates, we have stud- 
ied the behavior of the self-energy along the FS when 
approaching the DDW instability. It has been found 
that self-energy effects and spectral functions are very 
isotropic along the FS even for doping close to 5c- 

In Sec. V we have discussed the nature of the exci- 
tations, which interacting with the charge carriers pro- 
duce the self-energy renormalizations. Charge excita- 
tions, dominated by collective effects, are the main con- 
tribution to a'^F{uo). As collective charge fluctuations 
live on a large energy scale, they are the responsible for 
the large self-energy effects at large energy, producing 
the reduction of the quasiparticle weight and transfer- 
ring spectral weight to the incoherent spectra at large 
binding energy. 

PIH method seems to be a suitable alternative for cal- 
culating spectral functions in the t — J model, moreover 
it can be used independently or as a complement to other 
calculations as well. 
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